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Abstract 

Boolean networks are an important class of computational mod¬ 
els for molecular interaction networks. Boolean canalization, a type of 
hierarchical clustering of the inputs of a Boolean function, has been ex¬ 
tensively studied in the context of network modeling where each layer 
of canalization adds a degree of stability in the dynamics of the net¬ 
work. Recently, dynamic network control approaches have been used 
for the design of new therapeutic interventions and for other applica¬ 
tions such as stem cell reprogramming. This work studies the role of 
canalization in the control of Boolean molecular networks. It provides 
a method for identifying the potential edges to control in the wiring 
diagram of a network for avoiding undesirable state transitions. The 
method is based on identifying appropriate input-output combinations 
on undesirable transitions that can be modified using the edges in the 
wiring diagram of the network. Moreover, a method for estimating the 
number of changed transitions in the state space of the system as a 
result of an edge deletion in the wiring diagram is presented. The con¬ 
trol methods of this paper were applied to a mutated cell-cycle model 
and to a p53-mdm2 model to identify potential control targets. 


1 Background 

A Gene Regulatory Network (GRN) is a representation of the intricate rela¬ 
tionships among genes, proteins, and other substances that are responsible 
for the expression levels of mRNA and proteins. The amount of these gene 
products and their temporal patterns characterize specific cell states or phe¬ 
notypes. Thus GRNs play a key role in the understanding of the various 
functions of cells and cellular components, and ultimately might help to de¬ 
sign intervention strategies for the control of biological systems. Recently, 
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practical applications in cancer systems biology such as the identification of 
new therapeutic targets has stimulated the development of computational 
tools that can help to identify new intervention targets. Experimentally, the 
interventions are realized by manipulating the wiring diagram of a system 
with the use of drugs or by gene knockouts to impact the dynamics of the 
system so that it is directed towards a desired state !0 l9||32[[33] , From 
the modeling perspective, the identification of intervention targets amounts 
to finding a set of relevant nodes and edges that can be used for performing 
interventions in silico. 

Many dynamic systems theory approaches have been used over the last 
decades to develop computational tools for analyzing the dynamics of GRNs. 
As a result, a large variety of models exists today. Boolean networks is a class 
of computational models in which genes can only be in one of two states: ON 
or OFF. BNs, and more general discrete models in which genes can take on 
more than two states, have been effectively used to model biological systems 
such as the yeast cell cycle network 1201, the Th regulatory network [22], the 


lac operon [30], the p53-mdm2 complex [l|[5j[25], A. thaliana [3], and for 
many other systems [2,6,10,11,27| 341. 

Boolean networks (BNs) as models for GRNs were introduced by S. Kauff- 
and R. Thomas 129 . BNs have been proposed as a framework 
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that does not rely on kinetic constants and therefore requires fewer param¬ 
eters to estimate, which simplifies analysis. Boolean canalizing rules were 
introduced by S. Kauffman and collaborators 14 and reflect the concept of 
canalization in evolutionary biology that Waddington pioneered in 1942 [31. 


Boolean canalization has been intensively studied from the network dynamic 

It has been shown that networks that use only 
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nested canalizing rules exhibit more stable dynamics compare to network 
using random rules [|15,23 . Furthermore, it has shown that each additional 


layer of canalization provides a degree of stability 0211 . The Boolean func¬ 


tions in published models tend to have many canalizing variables 14,23 


Boolean Networks 

A Boolean network can be defined as a dynamical system that is discrete 
in time as well as in variable states. More formally, consider a collection 
xi,...,x n of variables, each of which can take on values in the binary set 
{0,1}. A Boolean network in the variables x \,..., x n is a function 

F = (/!,...,/„) : {0,1}"—>{0,1}" 
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where each coordinate function /) is a Boolean function on a subset of 
{x\,... ,x n } which represents how the future value of the z-th variable de¬ 
pends on the present values of the variables. 

The dynamical properties of a Boolean network are given by the dif¬ 
ference equation x(t + 1) = F(x(t))\ that is, the dynamics is generated by 
iteration of F. More precisely, the dynamics of F is given by the state space 
graph S, defined as the graph with vertices in K n = {0, l} n which has an 
edge from x E {0, l} n to y E {0, l} n if and only if y = F(x). In this context, 
the problem of finding the states x E {0, l} n where the system will get sta¬ 
bilized is of particular importance. These special points of the state space 
are called attractors of a Boolean network and these may include steady 
states (fixed points), where F(x) = x, and cycles, where F r (x) = x for some 
integer number r > 1. Attractors in Boolean network modeling might rep¬ 
resent cell types |16|j or cellular states such as apoptosis, proliferation, or 
cell senescence 


12 


28 


Canalizing functions 

A Boolean function f(x \,..., x n ) : {0, l} n —> {0,1} is canalizing in the 
variable Xi with canalizing input value a and canalizing output value b if 
f(x i,..., Xi = a,... , x n ) = b. That is, once Xi gets its canalizing input, it 
by itself determines the output of the function regardless of the value of the 
other variables. The variable Xi is called a canalizing variable. 


Nested canalizing functions 

Let a be a permutation on the set {1,2,..., n}. The function f(x i,..., x n ) : 
{0, l} n —> {0,1} is a nested canalizing function (NCF) in the variable order 
a^i),..., Xpfy, with canalizing input values ai,...,a n E {0,1} and canaliz¬ 
ing output values b ±,..., b n E {0,1} if it can be represented in the form 


/ 


b\ if x a ^ = ai 

62 if x a(l) a li x a(2) = a 2 

^3 if x a( 1) a li x cr(2) 7^ a 2j x a(3) = a 3 


bn if x a( 1) f 1 a l i ■ ■ ■ i x cr{n) — a n 
k bn + 1 if x a(l) 7^ a l, • • • , x a(n) 7^ a n- 


Let a = {ai,..., a n } and /3 = {bi ,..., b n }. Then / is said to be {<r, a, /?} 
NCF if it is NCF in the variable order x a m ,..., x a ^ with canalizing input 
values a and canalizing output values j3. 
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Partially nested canalizing functions 


Any Boolean function f(x i, ..., x n ) : {0, l} n — > {0,1} can be represented in 
the form: 

b\ if x^i) = a\ 

b 2 if x CT (i) ^ ai,x 0 -( 2 ) = 02 


/ = 


bd if x a(l) ^15 • • ■ 5 x cr(d) 

„ 9 if 7^ a li ■ • • > x o(d) 7^ a d 


where either n = d, where is a terminal canalizing variable and g is a 
constant, or n < d where g{x a ( d+ i ),..., x CT ( n )) is a non-constant function and 
none of the variables x a ^ d+ i ^,..., x a ^ are canalizing for g. The integer d is 
called the nested canalizing depth of /. Such Boolean functions are called 
partially nested canalizing functions (PNCFs), see [l8j for more details. 


Layers of Canalization 


A Boolean function can be represented in different forms as a nested canaliz¬ 
ing function. A unique representation of the function is obtained by grouping 
the variables in layers of canalization 
uniquely written as 


21 . Every Boolean function can be 


f(x i,... ,x n ) — M!(M 2 (... ( M r _i(M r P c + 1) + 1).. •) + 1) + b, (1) 

where Mi = + a q)> Pc is a polynomial with no canalizing variables, 

and k = k\ + • • • + k r is the canalizing depth. Each variable Xi appears in 
exactly one of the M±, M 2 , ■ ■ ■, M r ,P c . The proof of this property is given 
in |9|. The number r in Equation [I] is called the layer number of /. 

Example 1.1. Consider the Boolean functions f\, f 2 , and with truth tables 
given at Table [lj The layers representation for fi is 

/i(xi, x 2 , x 3 ) = xix 3 (x 2 + 1) + 1, where Mi = xix 3 (x 2 + 1), 


for f 2 is 

f 2 {xi,x 2 ,x 3 ) = (xi+l)[x 2 (x 3 +l)+l]-|-l, where Mi = aq+1, M 2 = x 2 (x 3 +l), 
and for f 3 is 

f 3 (xi,x 2 ,x 3 ) = (xi + l)(x 2 + x 3 + l), where Mi = xi + 1, P c = x 2 + x 3 + l. 


Thus fi has layer number equal to 1, f 2 has layer number equal to 2, and f 3 
has layer number 1. The polynomial P c does not have canalizing variables. 
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Xl X 2 X 3 

h 

Xl x 2 x 3 

h 

Xl X 2 x 3 

h 

0 0 0 

1 

0 0 0 

0 

0 0 0 

1 

0 0 1 

1 

0 0 1 

0 

0 0 1 

0 

0 1 0 

1 

0 1 0 

1 

0 1 0 

0 

0 1 1 

1 

0 1 1 

0 

0 1 1 

1 

1 0 0 

1 

1 0 0 

1 

1 0 0 

0 

1 0 1 

0 

1 0 1 

1 

1 0 1 

0 

1 1 0 

1 

1 1 0 

1 

1 1 0 

0 

111 

1 

111 

1 

111 

0 


Table 1: Truth table for the Boolean functions in Example 0 


Definition of control actions 


This paper considers two types of control action: deletion of edges and con¬ 
stant expression of edges. An edge deletion represents the experimental 
intervention that prevents a regulation from happening. This action can be 
achieved by the use of therapeutic drugs that target a specific gene inter¬ 
action, see reference |5| where this type of control has been experimentally 
applied. A node deletion can be represented by the deletion of all of its 
outgoing edges. A constant expression or a constitutive activation of a node 
might result in aberrant cell proliferation and cancer, see [5] where the con¬ 
stant expression of cyclin G in the signaling pathway of p53 is reported as 
a signature of abnormal gene expression that leads to cancer. But constant 
expressions could also help to drive the system into a more desirable state, 
see 27 where constant expression of nodes have been proposed as poten¬ 


tial controls. As a proof of principle, this paper will consider the constant 
expression of an edge as a potential control action. 


Definition 1.2 (Edge Control). Consider the edge Xi —>• Xj in the wiring 
diagram W. For E IF 2 , the control of the edge Xi —> Xj consists of 
manipulating the input variable x ? ; for fj in the following way: 


Fj{x,Uij) = fj{x h ,..., (uij + 1 )x h .. 


For each value of u t j we have the following control settings: 

• When Uij = 0, U{j) = fj(xj 1 ,...,Xi,...,Xj rn ). That is, the 

control is not active. 


• When u it j = 1, Fj{x, Uij) = ,..., Xi = 0,..., x jm ). This is the 

case when the control is active and the action represents the removal 
of the edge x t —> Xj. 
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For simplicity, in Definition 1.2 we considered only edge deletions. To in¬ 
clude both the deletion and constant expression of an edge we could consider 
the following control function 


Fj(x, fi) = .., (u+. + u id + 1 )xi + ufj ,..., x jm ) where u ip ujj E F 2 . 

Then for each combination of ufj and uf J we have the following control 
settings: 

• For uf = 0 ,uf = 0, Fj{x, 0,0) = fj(xj l ,. .., Xj,. .., xj m ). That is, the 
control is not active. 


• For u t = 1 ,uf = 0, Tj{x, 1,0) = fj(x jx ,. ..,x t = 0,...,x jm ). This 
action represents the knock out of the node Xj. 

• For u~ = 0 ,uf = 1, ^(x,0,l) = fj (x 3l , ,Xi = l,...,x jm ). This 
action represents the constant expression of the node Xj. 

• For u~ = l,7i+ = 1, Tj{x, 1,1) = f j (x jl ,...,x i + l,...,x jm ). This 
action changes the variable Xi to its negative value and might not be 
a relevant case of control. 


2 Methods 

Eliminating state transitions through edge deletion and con¬ 
stant expression 

We avoid undesirable state transitions in the state space graph of a system 
of canalizing functions by means of edge deletion in the system’s wiring 
diagram. 

Let F = (/i,...,/ n ) : {0, l} n —> {0, l} n be a Boolean network and 
S = (V,, E s ) be the state space graph of F, where V s C {0, l} n is the vertex 
set of S and E s C {0, l} n x {0, l} n is its edge set. Suppose for u, v £ V s there 
is a directed edge {u, v} e E s which represents an undesirable transition. 
We eliminate the transition by deleting appropriate edges from the wiring 
diagram of the system, W = (V w , E w ), where V w = {xi,..., x n } and E w = 
V w x V w . The following is a sufficient condition for eliminating a transition 
from S through deleting an edge in E w . 

Method 2.1. Suppose xt E V w which takes input from x k E V w , i.e. 
{xfc,xt} E E w (we will also use the notation Xk xt). Let also Xk be a 
canalizing variable in f t , the functions that determines the state of xt in S. 
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If the following four conditions are met, then deleting the edge {xk, xt} from 
E w results in eliminating the transition {u, v} from E s : 

1. No variable in a more dominant layer assumes its canalizing input in 

u. 


2. Xk has canalizing input 0. 

3. The k-th entry of u is 1, i.e. [u]*, = 1. 

4- Xk has canalizing output that is the negation of the t-th entry of v, 
that is [v]t. 


The reason behind the first condition is that if any variable whose layer 
is more or equally dominant than Xk s layer has assumed its canalizing input 
in u, then replacing Xk with 0 will have no effect on ffs output. Deleting 
an edge has to impose change on the network for control and so the second 
requirement is needed since if [u]* = 0 already, then deleting the edge Xk —> 
xt will have no effect on the network. The third condition has a similar 
explanation. 

Similar sufficient conditions can be stated for eliminating a state space 
transition through constant expression of an edge, simply by replacing 0 
with 1 and 1 with 0. 

Node deletion can also be used for control through canalization. In that 
case, node deletion corresponds to deleting the outgoing edges from the 
deleted node and Method |2.1| can be applied to each one of them. 


Effect of edge deletion and constant expression on the state 
space 

We count the maximum number of state space transitions that can be 
changed as a result of deleting a single edge. 

Method 2.2. Let F = (/i,..., f n ) : {0, l} n —¥ {0, l} n be a Boolean network 
where ft is a PNCF of depth d in m variables in canalizing variable order 
1,2 ,,d. The deletion of the edge Xk —> xt results in 

(a) up to changes in the state space if k < d and Xk is in 

the r-th layer of ft, where are the numbers of variables in 

layers 1,,r, respectively; that is, the probability that any transition 
will be removed from the state space upon deletion of Xk —> xt is at 
most /2 n = (l)^i+^+-+4. 
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(b) up to 2 n ~ d ~ 1 changes in the state space if d < k < m, i.e. Xk is 
not canalizing; thus the probability that a particular transition will be 
removed from the state space upon deletion of Xk —> xt is at most 

To see how the bound is calculated, notice that when Xk —> xt is deleted, 
half of the transition table can potentially change (the other half had Xk = 0 
already). Of the remaining half, half contains the canalizing input of a 
variable in the most dominant layer and so Xk cannot cause change. Now 
half of the half only can possibly change but half of that has the canalizing 
input of another variable from a more or equally dominant layer to the one 
where Xk is, thus preventing Xk from causing change, etc. 

This upper bound remains the same when instead of deleting an edge, 
an edge is constantly expressed. 


General procedure for identifying control edges 


Below we provide a general procedure for identifying control edges based on 
the two methods we developed. Figure [I] further illustrates it. 

Given a Boolean network model for a biological system: 

1. Formulate a goal in terms of the part of the state space you wish to 
avoid, e.g. a fixed point or a cycle. If it contains more than one state, 
consider all states that are part of it. Choose to begin with one of 
them, v. 


2. For control via edge deletion, identify all l’s in the state u preceding 
v; for control via constant expression, identify all 0’s. 


3. Begin with, say, the leftmost 1 (or 0) in u. Observe the canalizing 
structure of the functions in your model. Then check if the conditions 
of Method 12.11 are satisfied. 


4. If some of the conditions of Method 2.1 are not met or you are looking 
for other control options, proceed to the next 1 (or 0) in u. 


5. If you wish to avoid a cycle or other trajectory that contains several 
states, you can repeat the above steps on all states in order to find all 
control options. 


6. If you find multiple edges as candidates for control, you may want to 
choose to delete the ones that have the smallest impact on the state 
space thus minimizing the side effects of edge manipulation. 
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Figure 1: Flowchart of the procedure for identifying edges as control candi¬ 
dates. 


Results 


We apply the control methods we developed to the Boolean models of two 
networks: a model of the human tumor suppressor gene p53 pathways 
and a mammalian cell-cycle network |8|. 
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p53-mdm2 model 

In [IT], a Boolean model, Eq. of the widely studied p53 pathway is 
built, where the external signal is dna.dsb, the DNA damage input. 


ATM n ext = Wipl(ATM + dna.dsb) 

P 53 next = Mdm2(ATM + Wipl) 

Wipl next = p53 

Mdm2 nex t = ATM(p53 + Wipl) 

The other variables are ATM, p53, Wipl, and Mdm2. When dna_dsb 
= 0, the state space has a single fixed point, (0000), corresponding to no 
stress. However, when dna_dsb = 1, i.e. the DNA damage input turns on, 
the state space contains a single cycle of length seven (Fig. [3]) and no fixed 
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points. The cycle represents cyclic variation in the expression patterns of 
all the four genes. We want to prevent this cycle from taking place through 
removing one or multiple transitions from it. 



Figure 2: Wiring diagram of the p53-mdm2 model from |17j. Pointed arrows 
correspond to positive regulatory interactions and blunt arrows to negative 
regulatory interactions. The edge deletion control candidate is given as a 
dotted edge and the two constant expression control candidates are denoted 
by thicker edges. 


(1000) -5“ (1100) -s- (1110) 

A 

Y 

( 0001 ) ( 0110 ) 

A 


( 0011 ) ■< - ( 0111 ) 



Figure 3: The cycle for the Boolean network of the p53 pathways in the 
presence of DNA damage. The state space is defined as (ATM p53 Wipl 
Mdm2). 

The wiring diagram of the model is presented on Fig. [2] By Method 
we identify that deleting edge p53—>Wipl (which also happens to correspond 
to deleting the node Wipl) has the effect of removing the following four 
transitions from the undesirable cycle in Fig. [3j 


2.1 
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( 1100 ) -> ( 1110 ), 
( 1110 ) -> ( 0110 ), 
( 0110 ) -> ( 0111 ), 
( 0111 ) -> ( 0011 ). 


The bold entries in the states correspond to the entries where conditions 
3 and 4 of Method 2T apply. As a result, the system has only a single steady 
state, (1100). Unfortunately, not all four proteins are inactive in it but we 
will see that this can be achieved through constant expression of the edge 
Wipl—»ATM, as it will be discussed later. 

We can also count the number of changes that deleting p53—)Wipl from 


the wiring diagram induces on the state space transitions: by Method 2.2 
there can be up to 2 4-1 = 8 changes and in fact we observe exactly as many, 
demonstrating that the bound from Method 2.2 is sharp: out of the 16 
states in the state space, eight contain p53 = 1. When applying the original 
update rules, the value of p53 in these eight states remains 1 (and so Wipl 
becomes 1), while after deleting p53—>Wipl, Wipl becomes 0 since now the 
output of its update rule is 0, thus causing a change in the state space. 

Constant expression of an edge is another strategy for removal of transi¬ 


tions in a state space graph. There are analogous conditions to Method 2.1 


for constant expression, obtained by simply replacing 0 with 1 and vice versa. 
For example, Wipl is a canalizing variable in the function of ATM with 
canalizing input 1 and canalizing output 0. Therefore, we can set the edge 
Wipl—»ATM to constant expression in order to remove the following tran¬ 
sitions from the undesirable cycle: (10001) —» (11001), (11001) -» (11101), 
and (00011) —> (10001). The result is a state space with fixed point (0000), 
corresponding to no stress as when dna.dsb = 0. Another option for control 
via constant expression is the edge Mdm2—>-p53 which also results in a single 
steady state, although this time it is (1000). 


Mutated cell-cycle network model 

As a second application, we consider Faure et al. [:8| who proposed a Boolean 
model of the cell cycle progression. We focus on the scenario when the tumor 
suppressor retinoblastoma protein Rb is absent as reflected in Eq. ©• The 
wiring diagram for that case is given in Fig. [4j Faure et al. J8| assume 
that the expression of CycD changes independently of the cell’s content 
and reflects the state of the growth factor. According to their model, the 
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mammalian cell cycle with a mutated phenotype will cycle through the eight 
states (Fig. [5J even when CycD is inactive. 


E2F = (CycA A CycB) V ( P 27 A CycB) 

CycE = E2F 

CycA = (E2F A Cdc20 A (Cdhl A Ubc)) 

V (CycA A Cdc20 A (Cdhl A Ubc)) 
p27 = (p27 A (CycE A CycA) A CycB A CycD) 

V (CycE A CycA A CycB) (3) 

Cdc20 = CycB 

Cdhl = (CycA A CycB) V Cdc20 V (p27 A CycB) 

Ubc = Cdhl V (Cdhl A Ubc 

A(Cdc20 V CycA V CycB)) 

CycB = Cdc20 A Cdhl 



Figure 4: Wiring diagram of the mutated (CycD and Rb inactive) mam¬ 
malian cell-cycle network. Pointed arrows correspond to positive regulatory 
interactions and blunt arrows to negative regulatory interactions |8j. Edge 
deletion control candidates are dotted and constant expression control can¬ 
didates are denoted by thicker edges. 

We propose four edges from the wiring diagram in Fig.[4]that can be used 
for control in order to avoid the cycle in Fig. [5j These edges were identified 
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( 10010110 ) 

A 


( 11010100 ) 


( 11110100 ) 


( 00001110 ) 

A 


Y 

( 11100100 ) 


Y 

( 00101011 ) -<■- ( 00100011 ) - ( 01100000 ) 


Figure 5: The cycle which persists in the absence of a growth factor in the 
mutated phenotype (Rb inactive). The state space is defined as (E2F CycE 
CycA p27 Cdc20 Cdhl Ubc CycB). 


following Method 2.1 applied on transitions in the cycle with the objective 
of eliminating the cycle and also leading the system to fixed point (s) where 
p27 = Cdhl = 1 as in a normal cell. The results are summarized in Ta¬ 
ble [2j Other attempts at control produced fixed points and/or cycles where 
p27 = Cdhl = 0. 


Control action 

Outcome 

Delete E2F—»CycE 

Single fixed point 
with p27 = Cdhl = 1 

Delete Cdhl—>-Ubc 

> 90% of states in component 
with fixed point where 
p27 = Cdhl = 1 

Express constantly 
CycB—>E2F 

Two fixed points, 

both with p27 = Cdhl = 1 

Express constantly 
Cdc20—>CycA 

> 90% of states in component 
with fixed point where 
p27 = Cdhl = 1 


Table 2: Summary of control results for the mutated cell-cycle network. 


Discussion 

In the Methods section, we noticed that the bound from Method |2.2| is sharp. 
This was demonstrated using the p53 model. In general, the exact bound 
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from part (a) of Method 2.2 is achieved when Ft(x, ut,k) = or %k (as a 


function), where Ft{x,ut,k) is the function obtained from f t by plugging in 
the canalizing input of the variables that are more or equally dominant to 

Xk- 

The bound from Method 2.2 can also help choose which edge to delete 
or constantly express when there is more than one option with the purpose 
of controlling the side effects resulting from an edge manipulation. If it is 
desirable to minimize the impact on the state space, thus avoiding possible 
negative side effects on the system, one should choose for control an edge 
whose input variable is in the least dominant layer possible in the target 
function. That is, if Xi and Xj are both canalizing variables in ft and Xj is in 
a less dominant layer than Xi, then one should choose to delete or constantly 
express the edge Xj —> xt since according to the bound of Method 2.2, the 
maximum impact of this control is smaller than if Xi —> Xt is manipulated. 

It is important to point out that Method 2.1 only guarantees that a 
certain transition will be avoided and one may be able to use this to remove 
a cycle from a state space as we did for the p53 model. However, the 
method does not guarantee that the system will not contain other cycles 
since removing transitions from a cycle destroys the cycle but may also 
create a different one; nor does it guarantee that the resulting fixed point 
will be exactly the desired one as it was observed in both applications. To 
find controllers that give the desired fixed points one could use the algebraic 
methods described in 
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This paper considers edge manipulations as potential control actions 
to avoid undesirable attractors. Control through edge manipulations in the 
wiring diagram of a network has been previously considered in |4|. Although 
the authors of |4| consider edge additions in the wiring diagram as control 
actions, that is, by adding new regulators to the existing set of regulators 
that help the system to transition into a desirable attractor. 


Conclusions 

The structure of the canalizing variables in a biologically relevant Boolean 
rule plays an important role in the control of Boolean networks. Special 
combinations of canalizing inputs can help identify network controllers and 
the canalizing structure of a Boolean function allows to estimate the num¬ 
ber of transitions that change after using the type of controllers proposed 
in this paper. Moreover, the hierarchy of the canalizing variables can be 
used for assessing the impact on the network dynamics as a result of a given 
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control. This paper exploits the canalizing properties of Boolean rules to 
derive a method that can be useful for identifying control targets for avoid¬ 
ing undesirable states. Additionally, it provides a method for assessing the 
impact of the controllers on the dynamics of the uncontrolled network. Thus 
these two complementary methods can help in the selection of appropriate 
controllers. Method 2.1 gives a practical way for identifying the potential 
edges to control in the wiring diagram of a network for avoiding undesirable 
state transitions. Method 2.2 on the other hand, provides a measure of the 
impact of an edge deletion onto the state space of a model and establishes 
that this impact differs significantly based on the canalizing properties of 
the nodes involved: an edge coming from a node with stronger canalization, 
represented in the model by a variable in a more dominant layer, has expo¬ 
nentially higher probability to change the state space than an edge from a 
node with weaker canalization, represented by a variable in a less dominant 
layer. Therefore, Method 2.2 is a useful tool for assessing the impact of the 
controllers identified by Method 2.1 on the dynamics of the system providing 
a way for selecting desirable controllers. 
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